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Various parameterizations for the orbits under local uni- 
tary transformations of three-qubit pure states are ana- 
lyzed. The interconvertibility, symmetry properties, param- 
eter ranges, calculability and behavior under measurement 
are looked at. It is shown that the entanglement monotones 
of any multipartite pure state uniquely determine the orbit 
of that state under local unitary transformations. It follows 
that there must be an entanglement monotone for three-qubit 
ure states which depends on the Kempe invariant defined in 
A form for such an entanglement monotone is proposed. 
A theorem is proved that significantly reduces the number of 
entanglement monotones that must be looked at to find the 
maximal probability of transforming one multipartite state to 
another. 
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I. INTRODUCTION 

Entanglement is at the heart of the studies of quantum 
computation and quantum information theory. It is what 
separates these studies from their classical counterparts. 
If we are to understand what new phenomena occur when 
we look at the true quantum mechanical description of 
nature as opposed to the approximations of classical me- 
chanics then we must understand how the quantum me- 
chanical description differs from the classical description. 
Entanglement is a measure of this difference. While en- 
tanglement between two parties is quite well understood 
B ill Hi I the entanglement within a quantum algo- 
rithm or in a state shared between many parties involves 
multipartite entanglement which is just beginning to be 
understood §1 ll). 

An integral part of the study of entanglement is de- 
termining the probability of transforming one pure state 
into another by Local Operations and Classical Commu- 
nication (LOCC). For two part systems this problem is 
solved, or at least reduced to the problem of finding the 
eigenvalues of a hermitian matrix, by Q |^ . For a NxM 
pure state the Schmidt decomposition tells us we can 
write 
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where the X] are in increasing order, J^i -^1 = 1: the |i) 
and are an orthonormal set of vectors in space A and 
B respectively, and n — min(A^, M). If we define 



then the highest attainable probability of transforming 
1^) to !</)), F(l^) ^ 10)), is given by § 



P(|V) ^ 1<^)) 
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The proof of this theorem is constructive so we can actu- 
ally write down the transformation that gives us |0) from 
I'f/'). For pure states of more than two parts no such nice 
theorem is known. The question of whether two three- 
qubit pure states can be transformed into each other with 
non-zero probability by LOCC has been solved by Diir 
et. al. Q but just getting a reasonable upper bound on 
that probability when it is a non-zero is unsolved. In this 
paper I attempt to make some progress towards solving 
this problem for three-qubit pure states and hopefully 
shed some light on how we might solve it for larger di- 
mensional spaces and more parts. 

One way to find P{\ip) \<j>)) is to look at the entan- 
glement monotones i?(|'0)) for the two states. For the 
duration of the paper "state" will refer to a pure state 
unless explicitly called a mixed state. An entanglement 
monotone, EM, is defined as a function that goes from 
states to positive real numbers and does not increase un- 
der LOCC. As a convention the value of any EM for a 
separable state is 0. For mixed and pure states of any di- 
mension and number of parts the following theorem holds 
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where the minimization is taken over the set of all EMs 
pO| . This can be seen by considering P(p — > p') as an EM 
for p. The problem is that this minimization is difficult 
to take since there is no known way to characterize all 
the entanglement monotones for multipartite states. We 
would like a "minimal set" of EMs similar to the Ek for 
the bipartite case in order to take the minimization. 

The situation for three or more parts is somewhat dif- 
ferent than for bipartite pure states. Firstly, generic 
M X M bipartite states have a stabilizer (i. e. the set of 
unitaries that takes a state to itself) of dimension M — 1 
isomorphic to C/(l)®^^~^ while pure states with more 
parts generically have a discrete stabilizer. States whose 
parts are not of the same dimension may have larger sta- 
bilizers but bipartite states are the only ones that always 
have a continuous stabilizer. Secondly, the generalized 
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Schmidt decomposition, however you choose to general- 
ize it has complex coefficients for pure states 
with 3 or more parts. This implies that generically these 
states are not local unitarily equivalent to their complex 
conjugate states (i. e. the state with each of its coefficients 
complex conjugated). Also, for bipartite pure states all 
the local unitary (LU) invariants can be calculated from 
the eigenvalues of the reduced density matrices but this 
does not hold for more parts. I will go into more detail 
about LU invariants in the next section. 

The structure of the paper is as follows, in section || the 
interconvertibility, behavior under measurement, symme- 
try properties, parameter ranges and calculability of two 
generalizations of the Schmidt decomposition of equa- 
tion (0) and the polynomial invariants (defined below) 
are looked at. In section III it is shown that the entan- 



glement monotones uniquely determine the orbit of mul- 
tipartite pure states and this is used to show that there 
must be an EM algebraically independent of the known 
EMs. A form for this EM is proposed and studied. Sec- 
tion IV discusses other monotones that must exist and 



their properties. Lastly, in section ^ a theorem is proved 
that significantly reduces the number of EMs that must 
be minimized over to get P{p p') of equation (Q). 



II. DECOMPOSITIONS AND INVARIANTS OF 
THREE-QUBIT PURE STATES 



Let 1-0) be a multipartite state in Tii ®Ti2 ■ ■ - ^Tin and 
let A^^ -.TLi-^TL'^ be Krauss operators for an operation 

on the hilbert space Hi with Efe4'^^^fc^ = ^^ and Z, 
is the identity acting on Hi- A (non-increasing) EM is a 
real valued function E such that 
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for any state lip), operation and space i where 



Pk 



(5) 



(6) 



This definition for pure states is taken from the definition 
for a general state in |ic| ]. One can always transform a 
state into product states and a product state cannot be 
transformed into anything but another product state so 
the value of an EM for a product state is chosen to be 
zero and all other states must have a non-negative value 
for the EM. Since can be a unitary operator or the 
inverse of that operator, equation (^) implies that all 
EMs must be invariant under LU. Hence, a first step to 
understanding the EMs is to look at the LU invariants 
that parameterize the set of orbits. 

There are many ways to find LU invariants for three- 
qubit states H jl| § fll some of which 
can be generalized to more parts and larger spaces but 
for now I will concentrate on the three-qubit case. The 



three sets of invariants I will look at in this section are 
the polynomial invariants [|13|, what I will call the diag- 



onalization decomposition |12 
maximization decomposition 



and what I will call the 



A. The Polynomial Invariants 

A general polynomial invariant Pcr,T ilip)) for a state of 
the form 



is written as 
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where a and r are permutations on n elements, repeated 
indices are summed and t stands for the complex con- 
jugate of t [|l^. If one applies a unitary to any of the 
qubits in and explicitly writes out Pcr,r (1^)) again it 
becomes apparent that P^.r dV')) is invariant. Of course, 
any polynomial in terms of the polynomial invariants 
Pg. (I'll))) is another polynomial invariant. In fact, it can 
be shown that all the polynomial invariants arc of this 
form. 

We know from |l^] that generic three-qubit states have 
a discrete stabilizer so the number of independent poly- 
nomial invariants is given by 

dim [C^ (g, (g, C'^] - 3 dim[S'[/(2)] - dim[C/(l)] -1 = 5 

(9) 

where the last —1 is due to the fact that we are using nor- 
malized states. The 5 independent continuous invariants 
are 
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where eoo — eii — 0, and eoi = — eio — 1 and again 
repeated indices are summed. /4 is the Kempe invariant 
referred to in the abstract. If one writes out /s and uses 
the identity eye,s = SirSjs — SisSjr it can be shown that 
is just the sum and difference of 64 polynomials of the 
form in equation (ph. With one more discrete invariant. 



sign[lm[P( 
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the LU orbit of a three-qubit state is determined uniquely 
1^ 0. In this paper I will define sign [a;] as 1 for non- 
negative numbers and —1 otherwise. The polynomial 
invariants have the advantage of being easy to compute 
for any state and the four previously known independent 
EMs are the following simple functions of Ii, I2, I3 
and I5 



T{AB)C = 2(1 - h) 
'T(AC)B = 2(1 - I2) 
T{BC)A = 2(1 - I3) 

TABC — 2-\//5- 



(12) 



B. The Diagonalization Decomposition 

The diagonalization decomposition, DD, introduced by 
Acin et. al. [0 is accomplished by first defining matrices 
iTo)j,k = ioj,fc and (ri)^,^ = ti,],k, then finding a unitary 
operation on space A that makes Tq singular, finding 
unitaries on space B and C that make Tq diagonal and 
using the remaining phase freedom to get rid of as many 
phases as possible. What is left is a state of the form 



IV-dd) 



/Ti^iooo) + 



/Mle^'^llOO) 
VMi|110> + 



/mI|iii> (13) 



where > 0, /io + /ii + ^2 + + f^i = 1 and < cj) < n. 
Note that generically there are two unitaries that will 
make To singular but it can be shown that only one will 
lead to 4> between and tt. If there is another solution, 
with (f) between tt and 27r exclusive, it is referred to as the 
dual state of IV'DD)- Some nice properties of DD are that 
there is a 1 to 1 correspondence with the orbits and there 
are a set of invertible functions between the parameters 
of the decomposition and the set of polynomial invariants 
given above. Namely, 

/i = 1 - 2//o(Ai2 +/i4) - 2 A 
/2 = 1 - 2fiaip3+fM) - 2 A 

^3 = 1- 2/io(Ai2 + /i3 + fM) 

h = l- 3[(^2 + /i3)(M0 - fM) + A*4(l - f^i) 

-/^2Ai3Mo + (1 - Mo)(A " AiiA/4)] 

h = sign[sin(0)^^^/ii^2Ai3/^4 

x(A-^4(l-2Aio + A'i)-M2M3)] (14) 

where A = /ii/i4 + ^2M3 ^ 2.^/11/12/^3/^4 cos(0) and if we 
define 

Ji^\ (l-h-h + h-l^/h 



J2 = ^[l-h+h-h- 2V/5 



J3 = - I + /1-/2-/3-2V/! 



J4 



J5 = 



1 / 5 
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-h - 27/. 
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then the coefficients are given by 



^4 



J4 + J5 ± \/T 
2(Ji + ,h) 

A, i^2,3,4 
Mo 

, ± J2 + -h ^ Ji 

1 - /io ± 

Mo 

n^n^ + ^^t^^3 - Ji 



2^/ n^fi^fi^fif 



sign [ (sin (0^)] 



le sign[y fi^ /ig [Ji - J2 J3 
-MJ2 + J3 + J4- (Mo )')]] (16) 



where T = ( J4 + Js)^ - 4(Ji + J4)(J2 + J4)(^3 + J4) > 0. 
The + and — solutions for the coefficients correspond to 
I^Dd) ^^'^ dual state. The inversion of the equations 
for Ii was done independently in ||l^. Note that their 
definition of I 4 is different from the one in this paper. 

Another nice property of the DD is that we can per- 
form an arbitrary measurement on it in space A and stay 
in the DD form. Since any measurement can be broken 
into a series of two outcome measurements [^8|, we can 
look at the two outcome measurement Ai and A2 where 
A\Ai + AIA2 = I. Using the singular value decomposi- 
tion we can write A.^ = UiDiV where V does not depend 
on i because the two positive hermitian operators 
and ^2^12 sum to the identity and therefore must be si- 
multaneously diagonalizable. The diagonal matrices, 
can be written as 











(17) 



where < x, y < 1 Since we are only concerned with 
what orbit the outcomes are in we may choose the Ui 
transformation. Also, matrices of the form 











(18) 



where ij^i and '02 are real numbers, commute with the Di 
matrices so the most general V can be written as 



a 



— Vl — o?e 



(19) 



where < a < 1 and Q is real. If we choose 
ya 
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and similarly for U2 with (x, y) replaced with 

— a;^, a/I — j/^) then in going from IV'Dd) 
^iI^DD) ^'^^ coefficients undergo the following 

transformations 



Mo 

1 

Ml ^ - 

7 
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i = 2,3,4 



arg 



e-*'^(a;2 - y^)^^^^^^ _ ^2) + e'^^^ (21) 



and again similarly for ^2|V'Dd)' Things become more 
complicated when becomes larger than tt and we have 
a dual solution. In this case we need to transform to the 
dual state which can be quite tedious. It should also be 
noted that if we want to plug the new form for the DD 
coefficients into equations (^4|) then the normalization 
must be taken into account. The normalization will just 
be the sum of the new forms for fig through /i4. 



C. The Maximization Decomposition 

The Maximization Decomposition MD, has a 

somewhat different way of decomposing the three qubit 
states. First we find the states, |<?!>a), |0_b) and |0c) each 
defined up to an overall phase, that maximize 



5(|</.A),|0i3),|</>c)) = ||(^|0A)|0i3)|</'c)|l 



(22) 



and apply a unitary such that \(j)A)\(j)B)\4'c) becomes 
|000). Defining up to an overall phase, as the vector 
perpendicular to |0), then the derivative of g along |1) at 
the point |000), 



lim 

£^0 



5(|0)+e|l),|0),|0))-.9(|0),|0),|0)) 



2Re[(V'|100)(000|7/')] 



(23) 



must be zero because g(|0), |0), |0)) is a maximum. Since 
we still have phase freedom in |0) and |1) this implies that 
(■01100) = and similarly for (0|O1O) and (01001). Using 
the remaining phase freedom in the choice of |0) and |1) 
we can eliminate all but one phase leaving us with 

IV'MD) = ae^'^lOOO) + 6|011) + c|101) + d|110) + /|111) 

(24) 



where a^+b'^ 



-(P+P = 1, < < 27r, < a, 6, c, d, / 



and b,c,dj < a. Note that 5(|0a), |0b), |0c)) = 
Unfortunately, the parameters as they are given above 
are not in 1 to 1 correspondence with the orbits. While 
the decomposition is generically unique, there are choices 
of the parameters within the given ranges that are not 
the result of the decomposition. For example, states with 



+ e, 



have 



and any choice of 



g ^(|0) + |1)), ^(|0) + |1)), -=(|0) + |1)) > 



V2 



V2 
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(25) 



for e < 0.014. Hence, these choices of the parameters are 
not a result of the decomposition. The true ranges of the 
parameters that would give a 1 to 1 correspondence with 
the orbits are as yet unknown. 

A nice property of the MD is that is it symmetric in 
particle exchange. Exchanging the particles is equivalent 
to exchanging b, c and d. This makes the permutation 
properties of the polynomial invariants easier to see when 
written in terms of the MD coefficients. They take the 
following form 

I,^l-2{{a^ + d^){b^ + c') + a'f) 
I, = l-2{{a^ + c^){b' + d') + a^f) 
h = l-2{ia' + b')ic' + d') + a\f) 
A - 1 - 3(a2(l - a^) - {b^c^ + 52^2 + c2d2)(i _ 2a^) 

-2b^c^d'^ - 2abcdf cos(0)) 
/g = a2|a/2 + 46cde*"^P 

h &\gTL[abcdf sin(0)(a2(l - 2a^){l - 2c? - f) 

-Ab^c^d^ - 2abcdf cos(0))]. (26) 

It is apparent from these equations that /i, I2 and 73 
are symmetric in permutations of particles AB, AC and 
BC respectively and I a, I5 and Iq are symmetric in any 
permutation of the particles. Unfortunately, the equa- 
tions in (|2^) are not as easy to invert as those in (pH). 
In fact, just calculating the MD coefficients for an ar- 
bitrary state is not an easy task, as it is in the case of 
the polynomial invariants and the DD coefficients, since 
determining the unitaries for the MD involves maximiz- 
ing over a 6 dimensional space with typically many local 
maxima. 

One more interesting fact about the MD is that 1 — 
is a non-increasing EM. We know this because in ||l^ it 
is shown that a function of the form 



EkA,kB,kci\ip)) = ^ max \\Ta 



'^cWr (27) 



where Tx is a fcjf-dimensional projector on system X = 
A,B,C, is an non-decreasing EM and £^i.i^i(|-0)) = a^. 
The EM 1 — a? can be shown to be independent of the 
T from equation ( p^ ) by looking at the gradient vectors 
of the T, 1 - and N ^ a? + b"^ + + d'^ + p at, for 
instance, the point a — 3, b,c,d,f— 1 and (f) — ^. Since 
the gradient vectors span a 6 dimensional space, 1 — 
cannot be written in terms of the r and N. The problem 
with using 1 — a? as an EM is that one needs to find 
the global maximum of a 6 dimensional space with many 
local maxima to calculate it. This is a difficult task for 
most states. 
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III. FIFTH INDEPENDENT EM 



(JABC = 3 - (Ji + /2 + l3)h 



(31) 



In section || it was shown that all EMs must be in- 
variant under LU and hence are determined by the orbit 
of the state. For three qubit states this means that EMs 
are a function of only the polynomial invariants, DD co- 
efficients or MD coefficients. In fact, this determination 
is unique. 

Theorem 1 The set of all EMs for any multipartite pure 
state, uniquely determine the orbit of the state. 

Proof. Suppose two states and in Hi (^H2 ■ ■ ■ 
Tin have the same values for the EMs but lie in different 
orbits. We know by using equation that 



(28) 



so 1^) can be transformed to \(j)) (and vice versa) by n- 
party LOCC, n-LOCC, with probability 1. Since EMs 
are non-increasing with any n-LOCC they must remain 
constant during the entire transformation from lip) to |0) 
(and vice versa). Also, we know that any EM between a 
system X = A, B, . . . and the rest of the systems thought 
of as one (e. g. between B and {ACD...)), I will call 
these EMs 2-EMs, is also an EM for multipartite states. 
This is because any n-LOCC on the multipartite state is 
also a 2-LOCC between X and the rest of the systems, 
since the 2-EM is non-increasing over 2-LOCC it must 
also be non- increasing over n-LOCC. In particular the 
sum of the lowest k eigenvectors of the reduced density 
matrices, 



(29) 



(i. e. the 2-EMs in equation (||)) must be EMs. So the 
E^iliP)) must remain unchanged and hence the spectrum 
of px is unchanged during the transformation from lip) 
to \4>). In particular a measurement on space X, given 
by ^1 and A2, must be such that 



px 



(30) 



where N is the normalization. The only way this can be 
satisfied is if is a unitary matrix. This means that 

IV') and \4>) are unitarily equivalent which contradicts our 
original supposition. □ 

Since we know there are 5 parameters that determine 
the orbit of a three qubit state then by theorem 1 there 
must be 5 independent, continuous EMs. To the best 
of the author's knowledge the only 4 known independent 
continuous EMs that don't require a difficult maximiza- 
tion over a multidimensional space are the four r EMs 
defined in equation (12). Any candidate for the fifth inde- 
pendent EM must depend on I4 since the t are invertible 
functions of /i, I2, I3 and I5 respectively. The following 
function fulfills that criterion 



and numerical results suggest that it is an EM. After gen- 
erating over 300,000 random states and applying a ran- 
dom operation to each of them the inequality in equation 
(^) was never violated by aABC- Also, note that aABC 
is symmetric in particle permutations as is tabc- For 
the rest of the paper I will assume that aABC is an EM. 
Indeed, it may be that there is a set of measure zero or 
perhaps just a very small measure for which aABC is not 
a monotone and my numerical test didn't explore this 
space but there must exist some function of the polyno- 
mial invariants which is independent of the rs and is an 
EM. For it to be useful in improving our upper bound 
for -P(IV') — > 10)) there should be pairs of states lip) and 
16) such that 



fTABc{\(t>)) 



rim 



(32) 



and I have found such states numerically. The largest 
value of 



(^ABcil'ip)) 

crABc{\4>)) 



^(1-^)) 



(33) 



that I found in my limited number of examples of was 
0.01 and I was able to find examples of states for 
which T (l^)) /t {\(f))) is greater than one for all t and 
'^ABc{H))/<^ABci\(t>)) is less than one. 



IV. OTHER EMS AND THE DISCRETE 
INVARIANT 

The five independent continuous EMs, T(^ab)Cj '''(ac)b, 
'''(BC)Ai Tabc and cfabc, can easily be inverted to find 
Ii - /s but to completely determine the orbit of a state 
we must also have an EM that will give us the value of 
the discrete invariant /g. This is equivalent to finding 
an EM that is not the same for a state and it complex 
conjugate state. Note that Ii, . . .I5 and hence the r and 
<^ABC do not change when a state is conjugated but by 
looking at any of the sets of LU invariants we can see that 
generically a state is not LU equivalent to its conjugate. 
By looking at equation (^) we can sec that this implies 
that there must be EMs that are not the same for the 
generic state and its conjugate. It is also easy to see that 
for any operation that takes a state {ip) to its conjugate 
IV') with probability p there is an operation that takes 
IV') to IV") with the same probability. So, for a generic 
state IV') there must be an EM that goes down for the 
operation IV') {ip) and a similar one that goes down 
the same amount for |V') \ip)- So, EMs of the following 
form must exist 



(1^)) = 



V + V 
V 



o.w. 



1 



(34) 
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where v and v' are functions of T(^ab)Cj T(AC)b-: T(^bc)A, 
TABC and aABC- 

Also, from we know that there are two classes 
of three-part entangled states (i. e. states with 
'''{AB)Ci'''{AC)Bi'''(BC)A > 0) that Can be converted into 
each other with some non-zero probability within the 
class and zero probability between the classes. Namely, 
the GHZ-class which contains 



1 



|GHZ) = -j= (1000) 



(35) 



and has non-zero tabc and the W-class which contains 
|W) --^(1001) + 1010) + I100)). (36) 

and has tabc = 0- Looking again at equation (|^) we 
see that tabc tells us that -P(lV'w) ~^ I^GHZ)) ~ ^ 
but none of the previously defined EMs tell us that 

^(IV'GHz) ^ iV'w)) = is to have an EM that is fi- 
nite for GHZ-class states and infinite for W-class states 
or zero for GHZ-class states and non-zero for W-class 
states such an EM must exist. 



V. FINDING A MINIMAL SET 

Since T(^ab)c-, T{AC)b, T(^bc)a, tabc, ctabc and de- 
termine the orbit of the state all other EMs must depend 
on them. A fairly general way to create new EMs from 
known EMs is to use what I will call /-type functions 



Definition 1 A function f : S <Z W 
function if it satisfies the following 

1. /(O) = 



is an f-type 



2. if Xi > Ui for all i — 1,2, .. .n then f[x) > f{y) for 
x,yeS 

3. fipx + (1 - p)y) > pf{x) + (1 - p)f{y) for any 
x,y e S and < p < 1 . 



For a set of EM, {-E,}, we have 
'Ail^) 



Erim > pE, 



VP 



^^-HM^- '''' 



for any measurement Ai, A2 and any state ji/')- So, we 
have 



imm > f 
>pf 



pE 
E 



Mm 

VP 
VP 



{l-p)E 



+ il-p)f 



A2W 

V^ 



E 



V^ 



(38) 



where the first inequality comes from property 2 and the 
second comes from property 3. Hence, f{Ei, . . . , Em) is 



also an EM. We can show that any EM f{Ei, . . . ,Em) 
that is an /-type function of monotones Ei, . . . Em does 
not modify the upper bound on P{\'ip) — > \4')) given by 



~^ \v)) ^ nrm ■ 



(39) 



E^{m' 

First for the one dimensional case. 

Lemma 1 If f{x) is an f-type function with n — 1 then 

m^^;J-A (40) 



/(.)-"^n?' 



for any x,y eS. 
Proof. 

Case 1 For x > y from property 2 we know f{x) > f{y) 
and hence 



f{y) - 



(41) 



Case 2 For a; < y if we choose p = ^ e [0, 1) then we 
know from properties 1 and 3 that f{py) > pf{y) 
and so 



/(y) y 



(42) 



For n dimensions we have the following theorem 
(proved with S. Daftuar and D. Whitehouse). 

Theorem 2 If f{x) is an f -type function then 



^^^^ >min|-i,l!> 1 = 1,2,. ..n (43) 



for x,yeS. 
Proof. Let 



(44) 



then we have 

Case 1 If c > 1 then from property 2 f{x) > f{y) and 



so 



fiv) - 



Case 2 If c < 1 then define 



Zi — — I = 1, 2, ... n 

c 



(45) 



(46) 



and g{r) = f(rz). Notice that g{r) is an /-type 
function with n = I and hence 
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9ic) 



> c 



5(1) 

or substituting in / we have 



m 



> c. 



Using Zi > Ui and property 2 we have 



fiy) 



> c. □ 



(47) 



(48) 



(49) 



For three- 
qubit states if we take the minimum of E{\il!))/E{\(l))) 
over £ = {t(ab)c, T(ac)b, '^(bc)/!, t^bc ctabc, f we 
are actually taking the minimum over the infinite set of 
all /-type functions of £. Although from theorem 1 we 
know that all EMs must be a function of £ it is possible 
that there exist EMs that are not /-type functions of £. 
These EMs could cause /'(IV') to be lower than 

the minimum of E{\'il;))/E{W)) over £. The EM men- 
tioned at the end of section is an example of such an 
EM. 



VI. CONCLUSIONS AND FURTHER RESEARCH 

Theorem 1 along with theorem 2 implies that there 
should be a (not necessarily finite) minimal set of EMs, 
M, for which all EMs for three-qubit states or similarly 
for any type of multipartite states are /-type functions 
of M. I conjecture that such a minimal set should be 
simple since the /-type functions seem to be a rather 
general way of creating EMs that are functions of other 
EMs. The difficult part seems to be finding the EMs that 
are minimal and showing that they are minimal. Using 
numerical results it seems that the r may be minimal. I 
looked at functions of the t that are almost but not quite 
/-type such as r^ '^^ and numerically tested whether they 
are EMs or not. None of them were EMs. I cannot say 
the same for aABC and definitely not for since I do 
not have an explicit form for the v. 

There is further research that may help these prob- 
lems. If one could invert the equations in (26) to write 
a, 6, c, d, f and (j) in terms of /i , . . . , /g that would allow us 



to calculate the EM 1 — a not to mention find the ranges 
for and calculate the values of a, 6, c, d, / and (p. The EM 
1 — could be used to replace aABC or perhaps as an 
addition to £ and may prove more useful than (tabc ■ As 
far as finding the minimal EMs and showing that they 
are minimal, the arbitrary measurement on the DD at 
the end of section ( IIB ) may be useful since it allows us 
to look at the value of /i , . . . , /e before and after an ar- 
bitrary measurement on an arbitrary state with far less 
parameters than if we didn't take out the LU freedom. 
Also, it may be able to tell us the maximal probabil- 
ity of transforming the general complex state \tp) to its 



conjugate state IVJ) and this is a crucial piece of infor- 
mation that is needed to calculate v' in equation (|34|). 
Unfortunately, most of these tasks involve trying to solve 
nontrivial equations or systems of equations with many 
variables which can be difficult or even impossible. 
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